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Abstract 

A theory of weakly-nonlinear low-temperature relaxational absorption of 
acoustic and electromagnetic waves in dielectric and metallic glasses is de- 
veloped. Basing upon the model of two-level tunneling systems we show that 
the nonlinear contribution to the absorption can be anomalously large. This 
is the case at low enough frequencies, u;ro(T) <C 1, where to{T) is the mini- 
mal relaxation time for two- level systems with the inter-level splitting ~ ksT. 
In dielectric glasses, the lowest-order nonlinear contribution is proportional 
to the wave's intensity. It is negative and exhibits anomalous frequency and 
temperature dependencies, AF/Fo oc [ujtq{T)]^^/^T^^ . In metallic glasses, 
the nonlinear contribution is also negative, and it is proportional to the square 
root of the wave's intensity and to the frequency. Numerical estimates show 

that the predicted nonlinear contribution can be measured experimentally. 
PACS number: 05.45. +b 
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I. INTRODUCTION 



Acoustic and dielectric properties of glasses were intensively studied during the last 
decades. The most remarkable resultS of these studies was the discovery of universal low- 
temperature properties that are only weakly dependent on the chemical composition of 
the glass. These "anomalous" specific properties include low-temperature specific heat, 
thermal conductivity, thermal expansion, propagation of ultrasound, dielectric loss, electric 
and acoustic echo and some other properties. An important step towards an understanding 
of the "anomalous" properties of glasses was the introduction of the tunneling modelii. 
According to this model, there exist two-level systems (TLS) associated with local tunneling 
states in double-well potentials. These states are characterized by an energy difference 2A 
between the minima and a tunnel splitting 2A, the energy spacing being 



2e = 2VA2 + A2. (1) 

The parameters A and A are random, the distribution of A and In A being assumed constant. 
There are several review articlesi^ where experimental data and their interpretation on the 
basis of the tunneling model are given. 

Two different absorption mechanisms are usually discussed in connection with the TLS. 
The first one is resonant absorption, which is a direct absorption of acoustic (or microwave) 
quanta hu accompanied by transitions between the levels of TLS. Resonant absorption of 
an acoustic wave can be expressed asi 

r^'-^'^) = a{uj/s) tanh {huj/2kBT) (2) 

where a is a dimensionless coupling constant, s is the sound velocity, while T is the temper- 
ature. The case of microwave differs from Eq. @ only by the coupling constant, so we will 
discuss below only the case of acoustic waves. At high enough temperature {hu -C ksT) 
Eq. (D yields 

r(=^"^) ^ ahuo^/2skBT . (3) 



Another mechanism - relaxational absorption - is due to modulation of the populations of 
the TLS levels by the alternating deformation field created by the acoustic wave. This 
modulation is due to the time-periodic variation in the inter-level spacing 2e. This variation 
produces, in its turn, a modulation of the levels' populations which lags in phase the variation 
of e. This lag leads to the energy dissipation. The linear m. the acoustic intensity relaxational 
absorption can be estimated a^'i 



where ro(T) is the minimum relaxation time of the TLS having e ~ fc^T. Comparing Eq. 



If uJTo > 1 the ratio r^-'^^'^Vr^''"'^ ~ hu^To/kBT can be either greater or less than one under 
experimentally accessible conditions. 

Nonlinear absorption was mostly studied for the resonant mechanism. Actually, the 
nonlinear resonant absorption of sound is one of the basic effects confirming the existence 
of TLS in disordered materials (see e. g. Ref. |^). The factor tanh{huj /2kBT) in Eq. 
has a transparent physical meaning - it is the difference in equilibrium populations of the 
lower and upper levels of a resonant TLS having the inter-level spacing hu. The reason for 
the nonlinearity is that intense sound equalizes the levels' populations, which again results 
in a decrease in the absorption. Experimentally, the effect takes place at very low acoustic 
intensitiesi. Consequently, the resonant absorption decreases with the intensity increase, 
and at large enough amplitudes only the relaxational contribution can be observed. As a 
result, the non/mear relaxational absorption becomes important. 

The source of the nonlinear relaxational absorption is a strong modulation of the inter- 
level spacing of the relevant TLS (with the inter- level distance ~ fc^T) by the sound wave. 
As a result, the relevant TLS are able to absorb energy only during a part of the sound 
period. Consequently, the total absorption decreases with the sound amplitude. A theory 
for the nonlinear relaxation absorption in insulation glasses has been developed in0'§, the 
authors were concentrated on the case of large acoustic intensities. It seems that it is difficult 





with Eq. (^ we conclude that the relaxational absorption always predominates at ujtq ^ L 
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to realize such a regime in a realistic experimental situation with insulating glasses. How- 
ever, the nonlinear relaxational absorption was observed experimentally in metallic glasses 
(PdSiCu, PdSi, PdNiP)lil as a specific two-stage nonlinear behavior (successive decrease of 
resonant and then relaxational absorption as the intensity increases). The relevant theory 
for the case of large intensities has been developed in Refs. O. 

In this paper we show that one can expect pronounced effects in nonlinear relaxational 
absorption even at relatively low intensities. Namely, we claim that at low enough acoustic 
frequencies the nonlinear corrections to the linear absorption coefficient Tq^^^ are rather 
pronounced and possess unusual frequency and temperature dependencies. We hope that 
the regime of "weak nonlinearity" is easier accessible for the experiments, especially with 
insulating glasses. 

According to the theory for the strongly nonlinear regime0'0, significant deviation from 
the linear relaxational absorption takes place at d ^ ksT where d is the amplitude of 
modulation of the inter-level spacing (~ /c_bT) of the relevant TLS. Consequently, one can 
think that small nonlinear corrections to the linear absorption coefficient behave as (rf/fc^T)^. 
However, the energy difference 2e between the levels depends on both the diagonal splitting 
A and the tunnel coupling A according to Eq. ([I|). Furthermore, the relaxation time r for a 
given TLS is also dependent on its parameters A and A. This dependence can be expressed 
as (cf with Refs. |,|) 

1 1 / e /A\2 1 / e \ 

- coth — - . (5) 



r(A,A) ToiT)\kBTj \ej KksT , 

The first factor has a meaning of a minimal relaxation time for the systems with the inter- 
level spacing e ~ ksT . The fact that the inter-level transitions can take place only if the 
wells are coupled is described by the factor (A/e)^. This factor is maximal at A = or at 
A = e. In the factor (e/fc^T)^"'"^ one power of e is due to deformational coupling between TLS 
and internal degrees of freedom^i, while describes the energy dependence of the density 
of states for the internal degrees of freedom. In insulating glasses these are phonons {(3 = 2) 
while in metallic glasses the relaxation is due to electronic excitations near the Fermi level 



(/3 = 0)1. The last factor represents the sum of the occupation numbers for the excitation 
responsible for the absorption and emission processes. Indeed, coth.{e/ kBT) = 2N{2e) + 1 
where N is the Planck function. Consequently, nonlinear effects in absorption are determined 
by the ratio between the amplitude of modulation of the diagonal splitting, d, and the 
characteristic energy, S, which is a combination of the parameters A and A. This combination 
is determined by the dimensionless product ujtq, and it appears that 6 ^ ksT at ujt ^ 1 
(see below). As a result, at low frequencies nonlinear effects must be pronounced. They 
must also show specific frequency dependence. Having in mind that at ujt <^ 1 Tq^^^ oc u 
and temperature-independent, one can expect that nonlinear effects can be detected also 
from anomalous temperature and frequency dependencies of absorption. This is the main 
message of the paper. 

Below we will sketch the basic theoretical approach to calculate both linear and weakly 
nonlinear absorption. Then we will concentrate on the nonlinear correction to the absorption 
coefficient for relatively low amplitudes, both for insulating and metallic glasses. 

II. THEORETICAL BASIS 

The Hamiltonian of a TLS in the external ac field can be written aM 

TC = {A + d cos ujt)az — Aa^c , (6) 

where ax{z) are Pauli matrices, A is the energy gap between the levels of the isolated potential 
wells, while A is the tunnel matrix element, d = jikUik iii the case of sound wave {jik is the 
deformational potential of the TLS, is the amplitude value of the deformation tensor). 
In the case of electromagnetic wave d = rjSo, where rj is the dipole moment of the TLS 
while So is the amplitude of the electric field. The corresponding change in the tunneling 
transmittance A we shall assume to be small and neg lectili. 

According to the TLS model&i we assume A and A to be independent random quantities, 
A and In A being uniformly distributed over a broad range of values in comparison with the 
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temperature. We shall also assume the deformational potential to be a random quantity, 
uncorrelated with A and A , the distribution of which has a maximum (a similar assumption 
can be introduced concerning t]). 

Following Refs. we consider the case of relatively low frequencies when the energy 

of the acoustic quantum hu is much less that the characteristic inter-level distance S of the 
TLS which make a contribution to the nonlinear absorption (we shall estimate 6 later). At 
huj <^ 6 one can employ the adiabatic approximation and neglect time derivatives of the 
external field while solving the Schrodinger equation for the TLS. In this approximation the 
TLS is characterized by the time-dependent spacing 2e(t), 



e{t) = ^(A + c/ cos ujty + A2 , (7) 

and the occupation numbers of the upper (n) and the lower (1 — n) levels. The non- 
equilibrium occupation numbers can be found from the balance equation 

dn _ n - ripit) . . 

dt ~ T{t) ^ ^ 

where r(t) and no(t) are given by the substitution e{t) e into the Eq. and into the 
expression 

no = [exp{2e/kBT) + 1]"' (9) 

for the equilibrium occupation number. According Eq. (|]), the dynamics of the levels' 
population is characterized by a time-dependent relaxation time. The solution of Eq. (Rf) 



for a strongly nonlinear regime {d ^ ^sT) has been analyzed in Refs. p!0| , p!T| for the case of 
insulating glasses and in Ref. [T^ for metallic glasses. Below we shall present a theory for the 
regime of weak nonlinearity which seems easier to realize. 



III. CALCULATIONS 



In the following we will use the subscript to indicate linear results. The power absorbed 
by a single TLS can be determined by the expressionS 
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The contributions of individual TLS must be added and such a summation can be performed 
in a conventional way using the distribution function of the random parameters A and A 
and replacing the deformational potential by its average value. The distribution function 
is usuallyi'i chosen as 

A^(A,A) = 2iV/A 

where is the the density of states per energy interval while 1/A = (i(lnA)/(iA describes a 
smooth distribution of In A. The coefficient 2 is introduced because we use the notation 2e 
for the energy interval. This gives for the total absorption 

/•cxD fco AK 

P = 2Nj^ dAj^ — p(A,A,rf). (11) 

To analyze the absorption one can use the exact periodic in time solution of Eq. (|^) to 
obtainEi 

e Jo Jo ksTA V / 
e dtdt'e{t)e{t-t') 



Jo Jo cosh'^[e{t-t>)/kBT] ' ^ ' 

Here O = 271/^7 is the period of the wave. We will first consider the case of dielectric glasses. 
In the linear approximation one has 



eo(t) = -(A/e)u;rfsinLjt, e = VA^TA^, 

while the relaxation rate is time- independent and given by Eq. (|^). From Eq. (|T2|) one 
obtains the well-known expression for the relaxational absorption^ 



p Nuj r°° r°° dAdA cut 

~ 2kBT Jo Jo A cosh\e/ keT) 1 + cuV^ ' ^ ^ 



Since r(A, A) oc A^^ the typical value of A in this integral is ksT at utq ^ 1 and ksT ^ojtq 
at euro ^1. It is the last regime that produces anomalous nonlinear effects. In the linear 
case at cjTo ^ 1 one obtains 



Po = {'K'^/lQ)Nud^ (14) 

As we shall see, the nonlinear corrections to the absorption are of greatest importance for 
low values of the inter- level spacing e <^ ksT . Consequently, one can still use the asymptotic 
expression for coth(e/fcBT) ^ fc^T/e in Eq. (^). As a result, in the case oi (3 = 2 (i. e. for 
insulating glasses) r = TQiksT / h)'^ , and it is time-independent. That leads to important 
simplifications in Eq. (|12D. Indeed, to get nonlinear corrections to the absorption one has 
to expand only the function 

A + (i cos uot 
i[t) = = cud sin cut 

^{A + dcosLuty + A^ 

in powers of d up to second order, and then perform the integrations over t, t', A and A. 
Such a calculation can be easily performed using any computer algebra package (we used 
"Maple"). The result can be expressed as 

(P - Po)/Po = -Q {d/ksTf [uJToiT)]-'/' . (15) 

Here q = (3/32)(4-\/2 — 1) 0.44 is the numerical factor. 

In the case of metallic glasses, the linear calculations are performed analogous to the 
dielectric case. Making use of the A dependence of the relaxation time that emphasizes 
small A we get Pq ~ O-SSNood"^ for the most interesting case ujt <^ 1. For the nonlinear 
corrections, however, we can no longer assume the transition time to be time- independent. 
Indeed, calculations show that the main contribution to the nonlinear absorption comes 
from the expansion of the exponential in Eq. ([121) . Expanding the power p absorbed by a 
single TLS up to the fourth order in the wave's amplitude and subtracting the quadratic 
contribution pq (which is responsible for the linear absorption coefficient) we obtain 

z/ (i/ sin 9 — 8)n sin^ 9 cos^ 9 
(4 + z/sin^)(l + z/sin^) ^ ' 

Here we denote z/ = {ujtq)~^, A = esin6', A = ecos6'. 

As can be seen, this expression diverges strongly for small e, so it is not possible to obtain 
an exact quantitative estimate from our simplified calculation. In order to make an order-of 
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magnitude-estimate, however, we may cut off the integral in the lower limit, at e = \d\. This 
assumption yields for cur <^ 1 

(P - Po)/Po = -Cmiooro)-' {\d\/kBT) (16) 

with Cm is of the order 1. We observe that the nonlinear contribution is proportional to \d\, 
i. e. to the square root of the intensity rather than to the intensity, as in usual cases. This 
dependence arises from the fact that nonlinear effects are especially important for the TLS 
with small inter-level splitting e ^ \d\. Indeed, these systems are strongly perturbed by the 
acoustic wave. However, the probability to find a TLS with small splitting is proportional 
to \d\. Furthermore, at small e the relaxation time the ro(e) increases as that leads to 
a decrease of the contribution of the low-energy TLS. Collecting of all the factors leads to 
the estimate ([T6|). Unfortunately, our calculation based on the expansion in powers of d/e 
cannot provide the numerical factor. However, it yields the proper functional dependencies. 
In the high-frequency limit, ujt ^ 1, the nonlinear contribution is oc [d/kBT)"^. 

IV. DISCUSSION 

Since Pq (X uj and independent of T, the nonlinear contribution to the sound absorption 
in the dielectric glasses appears oc ^/LJ, the temperature dependence (at fixed d) being 
ocT-2[ro(T)]-i/2ocT-i/2. 

Analyzing the relevant integrals, one comes to the conclusion that the characteristic 
inter-level spacing of the TLS responsible for the absorption is 5 ~ kBT(ujToy^'^ < ksT. 
This fact confirms that one can employ low-energy asymptotics in Eq. (^) and replace 
cosh"^(e/A;Br) 1 in Eq. (|12]). Inequality hcu <^ S appears met for realistic parameters of 
the glass. 

Now, let us give a rough estimate of the intensities needed to obtain a significant non- 
linear contribution for dielectric glasses. One can estimate tq as {h/kBT){Tc/T)'^ where the 
constant depends on the interaction constant between TLS and thermal phonons. 
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As for the most glasses ~ 20 K we get at T = 0.1 K tq 0.3 • 10"^ s. Thus at 
iojl'K = 50 kHz lotq fa 10~^. We assume the deformational potential 1 eV and express 
the amplitude of the deformation tensor uf^ . The relation between the acoustic intensity 
W and the deformation tensor can be expressed as uf^ fa W/ (ps^) (where s is the sound 
velocity, while p is the density of the glass) . 

Demanding that the nonlinear correction should be about 10% of the absorption we get 



for p = 5 g/cm^, s = 3 • 10^ cm/s. 

In the case of metallic glasses we get Tq ~ h/xT with x ~ 0.01. This gives a total 
temperature dependence of compared to T~^/^ in the dielectric case. The characteristic 
interlevel spacing in the case of metallic glasses is 5 = {d ksT lotoY^'^ <S ksT. 

Using the same values for 7, p and s for the case of metallic glasses we get 



To conclude, we have estimated nonlinear contributions to low-temperature absorption of 
acoustic waves and microwave in dielectric and metallic glasses. It is shown that at low 
enough frequencies, u!To{T) <^ 1, they are anomalously large and can be detected experimen- 
tally. In the high frequency regime, ujto{T) > 1 the nonlinear contributions arc proportional 
to the intensity and do not exhibit anomalous frequency and temperature dependencies. 




fa 0.03[r(K)/0.1]^/^ 



(17) 



W (W/cm^) fa psWo (O.lfcsT/7)^ fti 0.01 [T (K)/0.1]^ . 
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